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Recent experiments have shown a remarkable number of coUapse-and-revival oscillations of the 
matter-wave coherence of ultracold atoms in optical lattices [Will et ah, Nature 465, 197 (2010)]. 
Using a mean-field approximation to the Bose-Hubbard model, we show that the visibility of coUapse- 
and-revival interference patterns reveal number squeezing of the initial superfiuid state. To describe 
the dynamics, we use an effective Hamiltonian that incorporates the intrinsic two-body and induced 
three-body interactions, and we analyze in detail the resulting complex pattern of coUapse-and- 
revival frequencies generated by virtual transitions to higher bands, as a function of lattice param- 
eters and mean-atom number. Our work shows that a combined analysis of both the multiband, 
non-stationary dynamics in the final deep lattice, and the number-squeezing of the initial superfiuid 
state, explains important characteristics of optical lattice coUapse-and-revival physics. Finally, by 
treating the two- and three-body interaction strengths, and the coefficients describing the initial 
superposition of number states, as free parameters in a fit to the experimental data it should be 
possible to go beyond some of the limitations of our model and obtain insight into the breakdown 
of the mean-field theory for the initial state or the role of nonperturbative effects in the final state 
dynamics. 
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I. INTRODUCTION 

Optical lattices have become a powerful tool for study- 
ing the quantum many-body states of ultracold atoms 
[TJI2]. In 1998 Jaksch et al. [3] suggested using ultracold 
atoms in an optical lattice to simulate the Bose-Hubbard 
model |4j, and in 2002 Greiner et al. |5| demonstrated 
Mott insulator quantum phase transitions by imaging 
atoms released from a lattice after a period of free expan- 
sion, which reveals the atomic momentum distribution 
at the moment of release. Recently, the in-situ measure- 
ment of the position of atoms has become possible [BHH]j 
and techniques to measure higher-order correlation func- 
tions, thermodynamic properties such as pressure, tem- 
perature, and transport properties have been developed 
[51[TT| . In the Mott state pairs of atoms have been stud- 
ied spectroscopically p!2Hl4] . 

In parallel, techniques have been developed to create 
and exploit the non-stationary dynamics of many-body 
states and their interactions. In collapse and revival 
experiments |15H17j . a sudden increase in lattice depth 
projects an initial ground state into a non-equilibrium 
state which is held in the deep lattice for a variable pe- 
riod of time t. Coherence between confined atoms at the 
moment the lattice is switched off results in an interfer- 
ence pattern in the atom-density after a period of free ex- 
pansion of the atom cloud [18]. The interaction-induced 
non-stationary dynamics in the lattice cause the degree 
of coherence to oscillate in time, which in turn makes the 



visibility of imaged interference patterns, defined here as 
the ratio of atoms in interference peaks divided by the 
total number of atoms, to oscillate as a function of t. 

In Ref. [in] we showed that the crucial role of multi- 
band physics in collapse and revival dynamics can be el- 
egantly described by extending the Bose-Hubbard model 
to include effective three- and higher-body interactions. 
After developing a new technique greatly reducing de- 
phasing from inhomogeneities. Will et al. [17| recently 
observed tens of collapse-and-revival oscillations, which 
confirmed the predicted higher-body effects [19]. The 
experimental data had characteristics, however, not ex- 
plained by higher-body interactions alone, and Will et 
al. [17] suggested that the data reveals additional infor- 
mation about the initial many-body state. 

In this paper, we show that a combined analysis of both 
the number-squeezing of the initial superfiuid state and 
the non-stationary multiband dynamics in the final deep 
lattice explains important characteristics of collapse-and- 
revival physics. First, we show how the visibility of the 
interference patterns depends on the degree of number- 
squeezing of the initial superfiuid in the shallow lattice, 
providing a new tool for analyzing many-body states in 
optical lattices. Next, to describe the dynamics in the 
deep lattice, we use an effective Hamiltonian that incor- 
porates the intrinsic two-body and induced three-body 
interactions, and we analyze in detail the resulting com- 
plex pattern of collapse-and-revival frequencies generated 
by virtual transitions to higher bands. In particular, we 
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show that the frequencies appearing the time evolution 
of the interference pattern visibiUties are robust under 
changes in mean atom number, with number squeezing 
only affecting the spectral weights. Finally, by treating 
the two- and three-body interaction strengths, and the 
coefficients describing the initial superposition of num- 
ber states, as free parameters in a fit to the experimental 
data it should be possible to go beyond some of the lim- 
itations of our model and obtain insight into the break- 
down of the mean-field theory for the initial state or the 
role of nonperturbative effects in the final state dynamic 

II. SHALLOW LATTICE AND INITIAL STATE 

The initial state of atoms in an optical lattice is well de- 
scribed by the ground state of the three-dimensional (3D) 
Bose-Hubbard model with its parameters in a regime 
where tunneling cannot be neglected. The many-body 
ground state to good approximation factorizes into a 
product of single-site states [5]. The breakdown of the 
factorization for small lattices has been studied in [20] . 

In previous analyses of collapse-and-revival experi- 
ments [TMT7] , the single-site state was modeled as a co- 
herent atom-number state = e^''*"'^ ''|0) with mean 
atom number (n) = (b^b) = where the operators b 
and b^ annihilate and create atoms in a lattice site, re- 
spectively, and |0) is the zero atom vacuum. This model 
captures basic aspects of the physics, but predicts con- 
trasts and spectral weights of the interference visibilities 
that differ significantly from the data (e.g. compare the 
theoretical predictions in [T^] to the data in TT]). For 
example, the data in |17j show minimal visibilities that 
are no smaller than 40% for some system parameters, be- 
fore dephasing from inhomogeneities is significant, which 
is much larger than the minimum visibility e~*^"^ pre- 
dicted for collapse of a coherent state. 

Below, we show that these deviations are associated 
with the number squeezing predicted by a mean-field 
decoupling approximation |21[ 122) . This approximation 
gives the Hamiltonian for each lattice site 

H, = -zJ,4){b + b^) + zJ,(t)'^ + hj2,,b^b^bh - ^ib'^b , (1) 

where z = 6 is the number of nearest-neighbor sites in the 
3D lattice, Ji is the tunneling parameter, U2,i is the on- 
site atom-atom interaction strength, and /i is the chem- 
ical potential, which sets the mean atom number. The 
ground state is the lowest eigenstate of the non-number 
conserving Hi corresponding to a real- valued <j) that mini- 
mizes the energy. For this value of (j) the superfluid order 
parameter (6) = 4>. A zero value for implies a Mott 
state with definite atom number per lattice site. Refer- 
ence |23j describes the breakdown of this approximation 
for ID and 2D systems. 

Explicit expressions for Ji and U2^i (see [22]) contain 
overlap integrals of Wannier functions with the single- 
atom Hamiltonian and the atom-atom interaction poten- 
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FIG. 1: (Color online) The probabilities |cnj^ for atom- 
number states |n) at each lattice site calculated from a mean- 
field decoupling approximation to the Bose-Hubbard model 
for several ^ — U2,i/{zJi) and fixed (n) = 2.46. The curve la- 
beled "coherent" corresponds to a coherent state with mean 
atom number (n). The inset shows the square of the order 
parameter as a function of ^. The values of ^ corresponds 
to a linear increase of the lattice depth from SEu to llEn 
assuming *''Rb atoms and laser wavelength of 738 nm. Here, 
En is the single-photon recoil energy of a ®^Rb atom. 



tial, respectively. We consider lattices where the Wannier 
functions can be approximated by harmonic oscillator 
solutions, and assume a regularized 3D delta function 
potential for the two-body interactions between atoms 
[HI with a strength proportional to the scattering length 
of atoms colliding with zero relative kinetic energy |25j . 
These approximations allow us to obtain useful analytic 
results, although fully quantitative modeling should in- 
clude the important role of the lattice anharmonicity. 

To obtain the ground state for different initial lattice 
depths we numerically solve for the ground state of Hi us- 
ing the variational solution l^f^) = J^n'^nl'n), with Fock 
states \n) containing n = 0, 1, 2, • • • atoms and ampli- 
tudes c„. We included n up to 20. Figure [l] shows exam- 
ples of probabilities |c„P for Fock state components |n) 
of the ground state for several ratios of ^ = U2^i/{zJi), 
but fixed (n). The probabilities for a coherent state \f3) 
with the same mean atom number are also shown. The 
ground state is number squeezed with increasing ^, and 
for large ^ the Fock states closest to (n) dominate. In 
the opposite limit ^ — it approaches a coherent state, 
however, number-squeezing is readily apparent even for 
a shallow lattice as long as ^ is large enough that a single 
well-defined band exists. The inset shows the superfluid 
order parame ter v ersus ^. For small ^ the parameter 
approaches while for large ^ the order parame- 

ter slowly approaches zero. We next analyze how the 
frequencies and visibilities of the collapse and revival os- 
cillations depend on the amplitudes c„. 
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III. EFFECTIVE DEEP LATTICE DYNAMICS 

Collapse and revival dynamics start after a sudden in- 
crease in lattice depth, reducing the tunneling param- 
eter to near zero. This operation is fast enough that 
interactions can be neglected but slow enough to pre- 
vent vibrational excitation in a lattice site. The state 
immediately after the lattice transformation retains the 
form X^n'^"!'^)' although the single particle basis func- 
tions change form slightly. Even with minimal real 
single-particle excitation, however, multiband effects are 
critical due to collision-induced virtual excitations to 
vibrationally-excited levels. We have shown that these 
effects can be concisely included via the effective Hamil- 
tonian [19] 

Hf ^ \^2ja}a}aa + ^-^jo} aaa O (Uij) , (2) 

where a^{a) creates (annihilates) an atom in a single 
renormalized spatial mode in a site of the final optical 
lattice and U2J and U^j give the strength of effective 
two- and three-body interactions. 

The effective Hamiltonian Hf applies to deep lattices 
with no tunneling, such that the dynamics in each site is 
independent. This is valid for the experiments in ll5H17j . 
where the tunneling and three-body recombination rates 
are 1-2 orders of magnitude smaller than the three-body 
interaction strength. Recently, Ref. [35] showed ways to 
measure tunneling rates from collapse and revival exper- 
iments in less deep lattices. Effective three-body inter- 
actions have also been shown to govern thermalization 
of a one-dimensional Bose gas |27| , and are important in 
Efimov physics [Hj. 

Renormalization fixes the value of U2J to the mea- 
sured energy of two atoms held in a lattice site after the 
zero-point energy of two non-interacting atoms is sub- 
tracted (Alternately, Ref. [21] gives the exact analytical 
relationship between U2J and the coupling constant of 
the regularized 3D delta function potential for two har- 
monically trapped atoms.) After renormalization, the ef- 
fective three- and higher-body interactions are predicted 
in terms of U2.f- Ignoring anharmonicities and assum- 
ing a harmonic potential with frequency ujf, we showed 
in Ref. [T^ that to second-order U^j is attractive and 
given by 

U3j = -CfUij/{h^f), (3) 

with Cf = 4v^ + 61n(4/(2 + ^3)) - 6 = 1.344.... The 
corrections due to anharmonicities in the deep lattices 
considered here are less than 5%. 

The state evolves under Hf for hold time t. Since 
Hf is number conserving, the state becomes |^'/(t)) = 
EnC„e-*^"*/'"», where 

En = U2jn{n - l)/2 + C/3,/n(n - l)(n - 2)/6 (4) 

and h — 2T:h is Planck's constant. 
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FIG. 2: (Color online) Visibility \{a)\'^ /\{a)t =o|^ as a function 
of hold time t after a sudden increase of lattice depth. The 
three panels correspond to three initial states characterized 
by the ratio ^ = U2,i/{zJi) before the increase. The mean 
atom number is (n) — 2.46 and the (renormalized) pair-wise 
interaction strength U2J = 0.0928fttj/ for all traces. Time 
is in units of h/U2j (e.g., 0.2 ms for *^Rb in a 41 Er deep 
lattice.) 



IV. INTERFERENCE PATTERNS 

After hold time t, the optical lattice is turned off and 
the atoms freely expand. The visibility of the spatial in- 
terference pattern versus t yields information about both 
the initial state and the dynamics. In [15 it was 
shown that the visibility is |(a)p/|(a)t=oP = I(q^)P/</'^j 
where the expectation value is over \^ fit)). A zero vis- 
ibility corresponds to a Gaussian spatial density profile, 
while a value of one corresponds to a diffraction pattern 
with maximum contrast. 

Figure [2] shows time traces of the visibility versus 
t for three initial states differing in the value of ^ = 
U2,i/{zJi), and (n) = 2.46 as in Fig. [1] The value 
U2.f = 0.0928hujf is typical for experiments with ^^Rb 
atoms in "deep" optical lattices. The traces show fast os- 
cillations with timescale ^ h/U2j, and a slower envelope 
with a timescale ~ h/U^j. The pattern is more complex 
for initial states close to an ideal coherent superfluid, and 
there are near total collapses of the matter-wave coher- 
ence at regular intervals. For increasing ^ the initial state 
is increasingly number squeezed, and the traces become 
smoother. In qualitative agreement with the data in |17) . 
there are only partial collapses of the visibilities. 

The value for (n) is chosen to accentuate the effect of 
both two- and three-body interactions and corresponds 
to easily obtained experimental parameters. For smaller 
(n) we find that the contrast and three-body effects are 
reduced, while for larger (n) the effective Hamiltonian 
needs to be improved. 

The behavior of the visibility in Fig. [2] can be eluci- 
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FIG. 3: (Color online) Spectrum of visibility as function of 
positive frequency v for initial coherent state (top panel) and 
two initial states obtained using mean-field decoupling ap- 
proximation (bottom two panels). Time traces of the latter 
two spectra are shown in Fig. [2| We use (n) = 2.46 and 
U2,f = 0.0928fta;/. Spectral lines correspond to the verti- 
cal lines. The spacings between dominant lines are shown in 
middle panel. The f = line strength increases from 0.2 to 
0.3 from the top to bottom panel. The arrow in the bottom 
panel shows the vanishing of a line at hv = U2J for large 
ratio U2,i/{zJi). For coherent states, the spectrum with and 
without effective three-body interaction are shown as full and 
dashed lines, respectively. For U3J = only frequencies with 
integer multiples of U2,f/h have nonzero line strength. The 
dotted curves are analytical expressions for the line strength, 
discussed in the text. 



dated by performing a Fourier expansion. The allowed 
frequencies v are differences between frequencies En/h 
governing the evolution of I ^'/(t)) [T7]. Spectra of the vis- 
ibility are shown in Fig. [3] The top panel shows the spec- 
trum for an initial, ideal coherent state, while the bottom 
two panels show spectra for states calculated with the de- 
coupling approximation. Consistent with the smoother 
traces in Fig. [2j fewer features are visible when ^ in- 
creases. For example, the arrow in the bottom panel of 
Fig. I3] highlights the disappearance of the line at U2j/h. 
This Dchavior is not obvious from examination oi Hf 
alone and dramatically illustrates why the analysis pre- 
sented here of the dynamics plus initial state is required 
to fully understand collapse and revival spectra. 

The frequency components of |(a)p follow a regular 
pattern, yielding a method for obtaining the two- and 
three-body coupling strengths from the dynamics. A pro- 
gression of equally spaced frequencies ends on or near 
each of the limits kU2.f/h for positive integer k. For 
example, for ly < U2,f/h features appear at 



U2j/h + mU3,f/h, 



(5) 



for m = 0, 1, 2, • • • . In fact, for v < kU2j/h and A: > 1 
the progression is kU2j/h + {mk + k{k — \)/2)Uy,j /h for 



m = 0, 1, 2, • • • corresponding to spacings between lines 
of kUsj/h. There are no lines at frequencies kU2j/h 
except for fc = 1. 

The Fock state amplitudes c„ of the initial state 
can, in turn, be directly obtained from the line strengths. 
For example, the lines at C/2,/ /h + mU^j /h have a 

strength of v^l™ + 2)(m -|- l)|c;,_^2Cm+id/'/'^ (the ex- 
pansion was previously shown in Ref. |17|). which de- 
pends on amplitudes of three Fock states. Again illus- 
trating that complete analysis of the spectra is subtle, we 
note that for m = the strength of the U2j/h spectral 
line is proportional to cq, the amplitude in the zero atom 
state! For example, for coherent or small ^ initial states 
the zero atom component has significant population (see 
Fig. fl]) and the line at U2j/h is visible; while for large 
^ it disappears. The variations in line strength there- 
fore correspond to the degree of superfluidity or number 
squeezing in the initial lattice state. 

For coherent states, closed- form expressions for the line 
strengths may be obtained allowing a detailed compari- 
son between the coUapse-and-revival spectra for coherent 
and more realistic squeezed states. For a coherent state 
and U3J = 0, the visibility is e2"^'('=°''('^2./*/'*)~i) [17], 
which in the frequency domain implies lines at kU2j/h 
for integer k with strength e^^*^ 7^(20^), where Ik{z) is a 
modified Bessel function. This spectrum is shown in the 
top panel of Fig. [3j where the dotted curve connecting 
the blue vertical lines is found by replacing the integer 
order k of the Bessel function by hv/U2j- 

The spectrum for an initial coherent state with U^j ^ 
can also be solved analytically. Frequencies U2j/h + 
mUsj/h have line strengths 



^2g-20^^4m/(p(^ + 2)r(m + 1)), 



(6) 



and frequencies 2U2j/h + {2m + l)U3j/h have line 
strengths 



(b*e-^'f''(j)^"'/{T{m + 3)r(m + 1)), 



(7) 



where T{z) is the Gamma function. Expressions for fre- 
quencies below kU2j/h for fc > 2 can also be found. This 
spectrum is shown in the top panel of Fig. [3j where the 
dotted curve connecting the lines is found by replacing 
m by the associated function of v. For frequencies below 
U2/h only lines near U2J /h+{n)U3j /h survive and, sim- 
ilarly, only frequencies near 2U2.f/h + (2(n) + l)U3j/h 
survive. 

We now study the dependence of the collapse and re- 
vival experiments on initial mean atom number, or equiv- 
alently the chemical potential. Figure |4] shows the visibil- 
ity as a function of hold time for several initial values of 
mean atom number. The system parameters Ji and C/2,i 
are held fixed and chosen such that the initial superfluid 
state shows significant squeezing. In contrast with Fig. 
where the patterns change dramatically with ^, Fig. 
shows oscillation patterns that are robust to changes in 
mean atom number. This robustness is because the fre- 
quencies in the time-traces do not depend on the average 



5 




123456789 10 
t (units of hAJj ^ 

FIG. 4: (Color Online) Visibility as a function of hold time 
for initial states with a mean atom number (n) of 1.5, 2.0, 2.5, 
3.0, and 3.5 going from the top to bottom panel. We use ^ = 
U2,i/{zJi) = 3.50 and f/2,/ = 0.0928fttJ/. The initial number 
probability and time trace for (n) = 2.5 closely resembles 
that for (n) = 2.46 and ^ = 3.50 shown in Figs, [l] and [2] 
respectively. 



bility distribution p((n)) over (n) and calculate 



p{t)= / d{n)p{{n})\{a)\\ (8) 
^0 



which is time dependent through (a). The visibility be- 
comes p{t)/p{0). Figure [5] shows time traces of the visi- 
bility after averaging over a Gaussian distributed p{{n)) 
normalized to one for (n) > 0. The three curves corre- 
spond to Gaussian distributions centered at an average of 
2.46 atoms, and with a square root of the variance (stan- 
dard deviation) that is either zero, 10%, or 20% of the 
average value, respectively. As in Fig. [4] we immediately 
note that the revivals occur at the same times with only 
small changes in amplitude, demonstrating that collapse 
and revival experiments are robust under small varia- 
tions in mean-atom number. A consequence of this in- 
dependence is that in the frequency domain there is no 
broadening of the spectral lines due to uncertainties in 
the mean atom number after averaging over experimen- 
tal realizations. There is broadening due to other mech- 
anisms. For example, inhomogeneities leading to a vari- 
ation of U2J across the lattice, as well as thermal effects, 
will cause both broadening and dephasing. The successes 
of the experiments in Ref. [.IZj are made possible by the 
ability to reduce these effects sufficiently to allow obser- 
vation of tens of oscillations. 



atom number, at least in the regime where the effective 
Hamiltonian Hf is valid. For example, the fast oscil- 
lations related to h/U2.f are in phase and the contrast 
is nearly independent of (n). This behavior can be un- 
derstood from our analysis of the frequency components, 
which are independent of c„ and thus also mean-atom 
number. (For example, for u < 1/2./ /h spectral lines ap- 
pear at C/2,/ /h + mils J //i for m = 0, 1, 2, • • • .) 

In Fig. |4] the differences in the time evolution are only 
noticeable at half integer multiples of h/U^j, correspond- 
ing to i sa A.0h/U2j in the figure. At these times the con- 



tribution to the phase factor 



appearing in the 



evolution of |\l//(t)) due to the three-body interaction is 
— 1 for the n = 3, 7, 11, ... Fock states |n). 




t (units hAJ^ j) 



V. MEAN ATOM-NUMBER UNCERTAINTIES 

We can quantify possible dephasing effects due to ex- 
perimental uncertainties in preparing the superfluid with 
the same mean atom number (n) or equivalently the same 
chemical potential. An experimental sequence requires a 
new preparation of the superfluid for each hold time as 
the measurement of the interference pattern is destruc- 
tive. We model this uncertainty by introducing a proba- 



FIG. 5: (Color online) Visibility as a function of hold time 
for three different averages over an experimental uncertainty 
in initial mean-atom number (or equivalently chemical poten- 
tial) in the superfluid. We use f = U2,i/{zJi) = 3.50 and 
U2,f ~ 0.0928/iaj/. The full (black) curve, already shown in 
Fig. [2] assumes complete reproducibility at (n) = 2.46, while 
the dashed (red) curve use a Gaussian distribution with a 
mean atom number of 2.46 and a standard deviation equal to 
0.25 (or relative uncertainty of 10%). The dash-dotted (blue) 
curve is obtained for a standard deviation of 0.5 (or relative 
uncertainty of 20%). 
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VI. CONCLUSIONS 

We have shown that coUapse-and-revival dynamics of 
atoms in optical lattices are strongly sensitive to both 
the initial many-body ground state and effective higher- 
body dynamics when projected into a deep lattice. The 
spectral analysis of the time evolution gives detailed in- 
formation about effective interaction strengths as well as 
amplitudes of atom-number states superposed in the ini- 
tial state. For example, number squeezing is apparent in 
the bottom two panels of Fig. [3] where fewer frequency 
components are present than for the coherent state shown 
in the top panel. We also found that the frequencies fol- 
low regular patterns that are independent of the initial 
superposition of number states, and the line strengths 
give information about the initial state. 

Our analysis is based on various approximations. Fu- 
ture work should evaluate corrections to the decoupling 
approximation for small lattice depths. For example, 
Ref. [23] showed a marked effect of dimensionality on the 
Mott phase transition, and Ref. analyzed the role of 
tunneling in collapse and revival dynamics. Finite tem- 
perature effects may also be important. Wedding-cake 



structures [3] that occur with additional weak confine- 
ment can be included using a local-density approxima- 
tion. Three-body recombination should be included in 
the analysis when in the future the dynamics are stud- 
ied over a longer period of time than currently observed. 
The inclusion of the anharmonicity of the lattice-site po- 
tential, which reduces the level spacing and affects the 
relationship between the effective two- and three-body 
interaction, can be taken into account to construct an 
improved theoretical model. 

Despite these technical limitations, however, the com- 
bination of effects considered here captures important 
features of the collapse-and-revival dynamics. By treat- 
ing the interaction strengths U2J and U^j, and the ini- 
tial state coefficients c„, as free parameters in a fit to the 
experimental data it should be possible to both test our 
model, and correct for some of its limitations, for exam- 
ple, the neglect of the effects of anharmonicities. Incon- 
sistencies in the fit could then provide insight into the 
breakdown of the mean-field theory (including thermal 
effects) for the initial state or the role of nonperturbative 
effects in the final state dynamics. 
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